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Abstract 

In this paper we have used our recently developed ah initio Wannier-function-based methodology 
to perform extensive Hartree-Fock and correlated calculations on Li20 and LiCl to compute their 
Born effective charges. Results thus obtained are in very good agreement with the experiments. In 
particular, for the case of Li20, we resolve a controversy originating in the experiment of Osaka 
and Shindo [Solid State Commun. 51 (1984) 421] who had predicted the effective charge of Li 
ions to be in the range 0.58 — 0.61, a value much smaller compared to its nominal value of unity, 
thereby, suggesting that the bonding in the material could be partially covalent. We demonstrate 
that effective charge computed by Osaka and Shindo is the Szigeti charge, and once the Born charge 
is computed, it is in excellent agreement with our computed value. Mulliken population analysis of 
Li20 also confirms ionic nature of the bonding in the substance. 

PACS numbers: 77.22.-d, 71.10.-w, 71.15.-m 
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I INTRODUCTION 



I. INTRODUCTION 

Ab initio calculation of dielectric response properties of materials are routinely performed 
using methods based upon density- functional theory (DFT).- Most of these calculations 
presently are based upon the so-called 'modern theory of polarization' which is based upon 
a Berry-phase (BP) interpretation of macroscopic polarization of solids.-"^ The practical 
implementation of the BP approach within various versions of the DFT, such as the local- 
density approximation (LDA), is straightforward because of their mean-field nature. How- 
ever, it is desirable to go beyond the mean-field level, so as to compute the infiuence of 
electron correlation effects on polarization properties such as the Born charge of crystals. 
Recently, we have proposed an approach which allows computation of various polarization 
properties of crystalline insulators within a many-body framework.-"^ The approach utilizes 
a Wannier-function based real-space methodology, coupled with a finite-field approach, to 
perform correlated calculations using a Bethe-Goldstone-like many-body hierarchy. Its suc- 
cessful implementation was demonstrated by performing ab initio many-body calculations of 
Born charges,- and optical dielectric constants,- of various insulating crystals. We note that 
an ab initio wave-functional-based real-space correlated approach has recently been success- 
fully demonstrated by Hozoi et a/.(author?)-to compute the quasi-particle band structure 
of crystalline MgO. 

In this paper we apply our recently-developed approach,^ to perform ab initio correlated 
calculations of Born charges of crystalline Li20 and LiCl. Li20 is a technologically impor- 
tant material with possible applications in thermonuclear reactors,- as also in solid-state 
batteries.- Based upon their infrared reflectivity and Raman scattering based experiment 
on Li20, Osaka and Shindo^ reported the value of effective charge of Li ions to be in the 
range 0.58 — 0.61, a value much smaller compared to the nominal value of unity expected 
in an ionic material. Therefore, they speculated whether the nature of chemical bond in 
the substance is partly covalent.- In this work we present a careful investigation of this 
subject, and resolve the controversy by showing that the effective charge reported by Os- 
aka and Shindo^ was the Szigeti charge and not the Born charge. Once the Born charge 
is computed from the Szigeti charge, its experimental value 0.95 is in excellent agreement 
with our theoretical value of 0.91, and much closer to the nominal value. We also present a 
plot of the Wannier functions in the substance as also the Mulliken population analysis to 
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confirm the ionic bonding in Li20. As far as LiCl is concerned, earlier we computed its Born 
charge at the Hartree-Fock level using our Wannier function methodology,- and found its 
value to be significantly smaller as compared to the experiments. Therefore, in this work, 
we perform correlated calculations of the Born charge of LiCl, and, indeed, obtain much 
better agreement with the experiments. 

Remainder of this paper is organized as follows. In section [Til we briefiy describe the 
theoretical aspects of our Wannier-function-based methodology. Next in section [1111 we 
present and discuss the results of our calculations. Finally, in section [TVl we present our 
conclusions. 



II. METHODOLOGY 



The Born effective charge tensor, associated with the atoms of the i-th sublattice, 

is defined as^ 

(1) 



f)piel) 



E=0 



where Zi is the charge associated with the nuclei (or the core) of the sublattice, VL is the 
volume of the unit cell, e is magnitude of the electronic charge, and Pa''^ is the a-th Cartesian 
component of the electronic part of the macroscopic polarization induced as a result of 
the displacement of the sublattice in the /3-th Cartesian direction, Awj/j. For small Awj/j, 
one assumes = , with /\Fa computed as the change m the electronic 

polarization per unit cell due to the displacement Awj^. For a given lattice configuration 
denoted by parameter A, we demonstrated that Pa''\X), for a many-electron system, can be 



computed as the expectation valued 



PtH>^) = ^{<''\X^^'^\<''), (2) 

Qe (= — e) is the electronic charge, (-^ oo), represents the total number of unit cells in 
the crystal, Xa'''^ = Xlfcii is the a-th component of the many-electron position operator 
for the Ne electrons of the crystal, and I^E'o^'*) represents the correlated ground-state wave- 
function of the infinite solid, expressed in terms of Wannier functions. For one-electron mod- 
els such as Kohn-Sham theory, or the Hartree-Fock method, we showed^ that Eq. ([2|) leads 
to expressions consistent with the Berry-phase-based theory of King-Smith and Vanderbilt^, 
and Resta^. 
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II METHODOLOGY 



Calculating the expectation value of Eq. ([2]) can be tedious for a general correlated 
wave function. We avoid this by using the finite-field approach whereby the expectation 
values of various components of dipole operator can be obtained by performing both the 
HF and correlated calculations in presence of a small electric field in that direction, and 
then computing the first derivative of the computed total energy with respect to the applied 
electric field.-"^ This approach of computing expectation values derives its legitimacy from 
the generalized Hellman-Feynman theorem, and is called "finite-field method" in quantum 
chemistry literature.— Therefore, the present set of HF and correlated calculations are dis- 
tinct from our previous ones^'^'^'^'^ in that the present calculations have been performed 
in presence of an external electric field. For details related to our Wannier-function-based 
HF methodology, in the presence of external electric fields, we refer the reader to our earlier 
papers.-"^ 

Correlation contributions to the energy per unit cell Ecorr, even in the presence of external 
electric fields, were computed according to the incremental scheme of Stoll^ 



where e^, Aeij, Aeijk, ■ ■ ■ etc. are respectively the one-, two- and three-body, . . . correlation 
increments obtained by considering simultaneous virtual excitations from one, two, or three 
occupied Wannier functions, and i, j, k, . . . label the Wannier functions involved.— Because 
of the translational symmetry, one of these n Wannier functions is required to be in the 
reference unit cell, while the remaining ones can be anywhere else in the solid. In our previous 
studies performed on ionic systems, we demonstrated that the increments can be truncated 
to two-body effects (n = 2), with the Wannier functions not being farther than the third- 
nearest neighbors.—*^ We follow the same truncation scheme here as well, with the method 
of correlation calculation being the full-configuration-interaction (FCI) as in our previous 
works, where other technical details related to the approach are also described.-'^'^'^ 



E, 



corr 
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III. RESULTS AND DISCUSSION 

In this section, first we present the results of our calculations of Born charges of Li20, 
followed by those of LiCl. 

A. Li20 

The primitive cell of Li20 consists of a three-atom basis (two Li and one O), with each 
atom belonging to an fee lattice. In our calculations anion was placed at the origin of 
the primitive cell, while the cations were located at the positions (±a/4 ± a/4, ±a/4). This 
structure is usually referred to as the anti-fluorite structure, to which other related materials 
such as Na20 and K2O also belong. Based upon intuition, one would conclude that the 
bonding in the substance will be of ionic type, with the valence 2s electrons of the two 
lithium atoms of the formula unit being transferred to the oxygen. Thus the two Li atoms 
will be in cationic state while the oxygen will be in dianionic state. However, the argument for 
partial covalence in the system goes as follows. It is well-known that the free ions do not 
exist in nature. Therefore, in lithium oxide it will be stabilized only because of the crystal- 
field effects. Still, one could argue, that ion in solid state will be very diffuse, leading 
to partial covalency. In light of these arguments, it is interesting to examine the results of 
the experiments of Osaka and Shindo^ performed on the single crystals of Li20, where they 
measured the frequencies of fundamental optical and Raman active modes using infrared 
and Raman spectroscopy. Using the measured frequencies and several other parameters, the 
authors fitted their data both to rigid-ion and shell models of lattice dynamics, to obtain 
values for Li effective charges of 0.58 and 0.61,- respectively. Noting that the obtained 
effective charge is quite different from the nominal value of unity, they argued that the 
result suggests that Li20 has partially covalent character.- Next, we argue, however, that 
the effective charge reported by Osaka and Shindo is actually the Szigeti charge,-!^ and not 
the Born charge. The shell-model value of 0.61 of the effective charge of Li was computed 
by Osaka and Shindo^ using the expression for the fiuorite structure derived by Axe^ 

^. = [(-lo-4oWoY/^^, (4) 

where vloI^to are the longitudinal/transverse optical phonon frequencies, tq = |, where 
a is the lattice constant, fi is the effective mass of the three ion system, while is the 
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III RESULTS AND DISCUSSION 



high-frequency dielectric constant of the material. The expression in Eq. [4] is nothing but 
the Szigeti charge^^ of the ion from which the Born charge (Zt) can be computed using the 
relation^ 

- + (5) 



3 ^ 

If we use the values eoo = 2.68 and Zs(Li) = 0.61 obtained by Osaka and Shindo^ in 
Eq. El we obtain the experimental value of Born charge of Li ions in Li20 to be 0.95 which 
is very close to the nominal value of unity. From this, the experimental value of the Born 
charge of the O ions can be deduced to be 1.90. These values of the Born charges of Li and 
O ions imply, in the most unequivocal manner, an ionic structure for Li20 consistent with 
the intuition, without any possibility of covalency. 

Next, we report the results of our ab initio correlated calculations of Born charges of 
bulk Li20. These calculations were performed using the lobe-type Gaussian basis functions 
which, as demonstrated in the past, represent Cartesian Gaussian basis functions to a very 
good accuracy.— i^"^ As in our previous HP study of the ground state properties of Li20,— 
our calculations were initiated by a smaller basis set reported by Dovesi et alM consisting of 
(2s, Ip) functions for lithium, and (4s, 3p) functions for oxygen. However, we also augmented 
our basis set by a single d-tjpe exponent of 0.8, recommended also by Dovesi et al.,"^^ so as 
to check the influence of higher angular momentum functions on our results. Henceforth, 
we refer to these two basis sets as (s,p) and (s,p, d) basis sets, respectively. For the purpose 
of the calculations, fee geometry along with the experimentally extrapolated T = lattice 
constant value of a = 4.573A^ was assumed. For our finite-fleld approach to Born charge 
described earlier, electric fields of strength ±0.001 atomic units (a.u.), along with an oxygen 
lattice displacement Au = O.Ola, both in the x— direction, were used. This leads to the 
determination of Z*^ component of the Born charge corresponding to the oxygen lattice. 
During the correlated calculations, Wannier functions corresponding to Is functions of both 
Li and O were held frozen, while virtual excitations were carried out from the 2s and 2p 
orbitals of oxygen anion. At the HF level we verified that the symmetry relation Z*^ = 
Zyy = Z*^ holds, along with the fact that the off-diagonal elements of the Born charge 
tensor are negligible. Additionally, the optical sum rule Z*^(0) + 2Z*^(Li) = 0, was also 
verified. 

Our calculated results for the Born charge of oxygen in Li20 are presented in table [B It is 
obvious from the table that irrespective of the basis set used, most of the contribution to the 
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Table I: Correlated values of Born Charge of \j\2O obtained in our calculations. Column with 
heading HF refers to results obtained at the Hartree-Fock level. Heading "one-body" refers to results 
obtained after including the corrections due to "one-body" excitations from each Wannier function 
of the unit cell, to the HF value. Two-body (O) implies results include additional corrections 
due to simultaneous excitations from two distinct Wannier functions located on the anion in the 
reference unit cell. Two-body (INN) and two-body (2NN) correspond to two-body correlation 
effects involving 1st and 2nd, nearest neighbor Wannier functions, respectively. 



Basis Set 


Oxygen Born Charge 


Experiment 




HF 


One-body 


Two-Body (0) 


Two-body (INN) 


Two-body (2NN) 


1.90'' 




1.794 


1.783 


1.763 


1.781 


1.777 




{s,p, d) 


1.851 


1.840 


1.814 









As deduced from the Szigeti charge reported in reference.— 



Born charge, as expected, is obtained at the HF level, with relatively smaller contribution 
of electron-correlation effects. Correlation effects do diminish the value of the Born charge 
as compared to its HF value, however, rather insignificantly. Comparatively, speaking the 
inclusion of one type basis function on oxygen has greater influence on the Born charge of 
oxygen, and its value is increased as compared to the one obtained using the (s,p) basis set. 
With the {s,p,d) basis set, we found negligible influence on the Born charge from the two- 
body contributions outside of the reference unit cell. Thus our calculations, as per optical 
sum rule, predict a value of Born charge of Li to be close to 0.91, in very good agreement 
with the experimental results,- with the error being less than 5%. In order to obtain a 
pictorial view of electronic structure in the material, in Fig. [T] we plot one of the oxygen 2p 
valence Wannier functions, computed using the (s,p, d) basis set, along the (1, 1, 1) direction. 
To simulate the effect of one of the phonon modes, we displaced the oxygen atom of the 
primitive cell to the position (O.Ola, O.Ola, O.Ola), where a is the lattice constant, leaving 
the Li atoms undisturbed. From the figure it is obvious that the Wannier function is highly 
localized around the oxygen site, with small value in the interatomic region, implying the 
ionic character of the material. 

Another quantity which could help us determine whether the system is ionic or covalent 
is the Mulliken population associated with the atoms of the primitive cell. With the {s,p) 
basis set, the Li and O Mulliken populations were found to be —2.005 and 0.984, respectively. 
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III RESULTS AND DISCUSSION 



Figure 1: One of the 2p valence Wannier functions of oxygen plotted along the (1, 1, 1) direction. 
In order to simulate the effect of optical phonons, calculation was performed with oxygen shifted to 
the position (O.Ola, O.Ola, O.Ola), and the Li positions unchanged. Nodes in the Wannier function 
at positions (±a/4, ±a/4, ±a/4) are due to its orthogonality to the Li Is Wannier functions located 
there. 




.il I I I I I I I I I I I I I I I I I I I I 

-5 -4 -3 -2 -1 1 2 3 4 5 

r (a.u.) 

while with the {s,p,d) basis set the populations were —2.004 and 0.984. It is obvious that 
the computed Mulliken populations have converged with respect to the basis set, and are 
very close to the nominal ionicities of the two atoms. This further confirms the ionic picture 
of the bonding in this substance. 

B. LiCl 

These calculations were performed using the basis set proposed by Prencipe et al.— in 
their ab initio HF study of alkali halides. The experimental fee geometry with lattice 
constant of 5.07 A was assumed. The unit cell employed during the calculations consisted of 
anion placed at the origin, and the cation at the position (a/2, 0, 0). During the many-body 
calculations, only the valence Wannier functions composed of 3s, and 3p orbitals localized 
on CI, were correlated. In order to calculate the Born charge, the anion was displaced in the 
x-direction by an amount O.Ola. For the finite-field computation of the dipole expectation 
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Table II: Correlated values of Born Charge of LiCl obtained in our calculations. Rest of the 
information is same as in the caption of table [B 



Basis Set 


Chlorine Born Charge 


Experiment 




HF 


One-body 


Two-Body (0) 


Two-body (INN) 


1.231'' 


{s,p) 


1.019 


1.020 


1.037 


1.013 




{s,p, d) 


1.098 


1.128 


1.178 







"Computed from the experimental data reported in reference.— 

value, an external electric field of strength ±0.01 a.u., also along the x-axis, v^as applied. 

The results of our calculations are presented in table [Til With our (s,p) basis set the 
correlation effects v^ere included up to the level of two-body increments involving the Wannier 
functions in the nearest-neighboring cells. However, with the augmented {s,p,d) basis set, 
inclusion of two-body increments beyond the reference unit cell, did not make any significant 
difference to the results. Similar to the case of, Li20, we note that: (a) inclusion of one 
d-type function improves the results both at the HF and the correlated levels, and (b) the 
inclusion of electron-correlation effects generally improves the results as compared to the HF 
values. The final value of Born charge of LiCl of 1.178 obtained from these calculations, is 
again in very good agreement with the experimental value, with an error smaller than 5%. 



IV. CONCLUSIONS 



In conclusion, we performed Wannier-function-based ab initio correlated calculations on 
crystalline Li20 and LiCl to compute their Born effective charges, using a methodology re- 
cently developed in our group. Calculations were performed using polarizable basis sets, and 
the results obtained were in very good agreement with the experimental ones. Additionally, 
we also resolved an ambiguity associated with rather small value of experimentally measured 
effective charge of Li20, which had raised the possibility of partial covalence in its chemical 
structure,- by demonstrating that the measured charge was the Szigeti charge, and not the 
Born charge of the compound. The Born charge obtained from the Szigeti charge was found 
to be in good agreement with its nominal ionic value. By further computing the Mulliken 
charges and exploring the Wannier functions, we demonstrated that the results were con- 
sistent with a ionic picture of Li20. The computational effort involved in our calculations 
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IV CONCLUSIONS 



is significant heavier than in the calculations performed using the density-functional theory. 
This is because our approach is many-body in nature, requiring calculations beyond the 
mean-field. At present we are exploring the possibility of using our approach to compute 
the Born charges of ferroelectric, as well as covalent materials, and we will publish our results 
in future as and when they become available. 
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